Analytic proof of the Sutherland conjecture 



J. Dittrich^ 

Nuclear Physics Institute ASCR, CZ-250 68 Rez, Czech Repubhc 

and 

V.I. Inozemtsev 
BLTP JINR, 141980 Dubna, Moscow Region, Russia 



Abstract 

Using the integral representation of the inverse of the logarithmic derivative of the 
elliptic theta function, the spectrum of the Lax matrix for the ID system of particles 
interacting via inverse sinh-squared potential is shown to be given by the asymptotic 
Bethe ansatz in the thermodynamic limit. 
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The problem of the verification of the asymptotic Bethe ansatz method [1] still re- 
mains unsolved more than 30 years after its first presentation. The method consists in 
using only scattering data for the description of the integrable many-body systems in the 
thermodynamic limit. It is well known that if the system is integrable in the Yang-Baxter 
sense, the many-body scattering matrix is expressed via only two-particle phase shift but 
real structure of the wave functions might be rather complicated if the interaction is non- 
local. In particular, when the two-body potential is of the form sin~^(7ra;/L), where L is 
the size of the system, or 1/ sinh^(a;), the wave functions differ drastically from the linear 
combinations of the plane waves inherent for the Bethe ansatz. Despite the exact results 
in the thermodynamic limit available for sin~^(7ra;/L) case are in complete coincidence 
with the asymptotic Bethe ansatz, the reason for this coincidence is still quite mysterious 
and it cannot be used as an argument to validate the method for the case of 1/ sinh^ x 
pair potential which is more complicated from the mathematical viewpoint . 

Due to the lack of a general approach to the problem, any particular exact results 
confirming the asymptotic Bethe ansatz are of interest. Some years ago, Sutherland [2] 
proposed one example of very good numerical coincidence of the asymptotic results with 
exact ones. It concerns the densities of the distribution of the eigenvalues of the L matrix 
from the Lax pair [3] for the systems of particles interacting via 1/ sinh^x potential, 

Ljk = PjSjk + (1 - Sjk)iX coth.{xj - Xk), 

where pj — —id/dxj obey the canonical commutation relations = iSjk- The corre- 

sponding Hamilton operator reads 

2 j=i j<fc smh [Xj - Xk) 

Asymptotically, if xi <^ ... <^ xn, the particles have the momenta ki < ... < kN and the 
elements of the Lax matrix become c-numbers, 

{Las)jk = kj6jk + i\sgn{j - k). (1) 

The asymptotic Bethe ansatz gives the asymptotic momenta as solutions to the equations 

Lkj^2nIj + J2r{kj-ki), (2) 
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where T{k) is the two-body phase shift, L is the total size of a system and {Ij} are 
quantum numbers. In the classical limit A — > oo, for the ground state, (2) becomes an 
integral equation for large N and L [4], 

rA 



where 



2a— dx'^{x — x')p{x'), 

J— A 



'j{x) = ln(l + X ^), 



(3) 



p{x) is the density of the momentum distribution in the ground state, and a — L/N is 
average nearest-neighbor spacing (or the lattice constant). The kernel of this integral 
equation is symmetric and positive definite. Thus it has unique solution at given A [6]. 
The normalization condition 



f 

/ p{x')dx' 

J -A 
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(4) 



defines ^4 as a function of a. The distribution of the eigenvalues of the Lax matrix can be 
connected with the distribution of the momenta [2]. Indeed, it follows from (1) that the 
equation for the eigenvalues 



det(Las - Iz) 



N N 



iX) 



N 



s=l 



can be written as 



N 



kj — Us + iX 
ikj — ujs — iX 



N 



— ^ arctan 



A 



ikj — Us 



A discontinuous branch of arctan with values in [0, tt] is used here. In the thermodynamic 
limit, the eigenvalues {uj} are distributed with the density (7{uj): Na{u!)du! gives the 
number of {uj} in the interval {00,00 + duo). Hence 

— = Naiu). 
duo 

Differentiating the above relation with respect to uo and taking classical limit gives 



1 r 



dxp{x) 



a{x- + 1/4 



(5) 
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after rescaling of variables [2]. One can see that in the classical limit the density a{ijj) 
can be calculated via the solution of the integral equation of the asymptotic Bethe ansatz 
method (3). 

On the other hand, in this limit the particles take their equilibrium positions at Xj = ja 
in the ground state with kj — 0. The form of the Lax matrix becomes simpler, 



Ljk = (1 - Sjk) i A coth[a(j - k)], 



(6) 



and the distribution of its eigenvalues can be calculated directly since (6) is of the Toeplitz 
form. Its eigenvectors are the plane waves, and after imposing periodic boundary condition 
(i.e. regularization of the determinant) and taking thermodynamic limit N ^ oo one 
could introduce the variable = 2t:s/N defining the continuous distribution 

The result can be written upon rescaling uj —>■ 2\uj in the form [2] 

9[{4>I2) 



a;(0) = 



(7) 



2^1(0/2)' 

where 9i{x) is the standard theta function 

oo 

e^{x) = 2 ^(-l)"5("+V2)^ sin(2n + l)x 

n=0 

with the nome q — e~". The density of the eigenvalues is given now by the formula 



where the derivative is calculated through the relation 



dcu 
d4> 



"7r2 



1 

~4 

K-E 



^i(0/2) 

e,{<t>/2) 



K sn2(X0/7r)J ■ 

Thus one gets two representations for the density of the eigenvalues of the classical Lax 
matrix, one exact ( formulae (7-9)) and one obtained by using the asymptotic Bethe 
ansatz method (formulae (3-5)). If it is true, they should coincide. The main difficulty 
in verifying this fact is that there is no chance to find analytic solution to the integral 
equation (3). In [2], Sutherland found good coincidence of both expressions by solving 



this equation numerically with high accuracy. However, analytic solution of the problem 
has not been found. 

In what follows, we propose a construction which uses analytic properties of the elliptic 
functions and provides the desired proof. Let us introduce the notation Xr = 
Xi — '^smx for any complex x- Consider at first the problem of explicit construction of 
the inverse function such that 

0(u;(A))^A. 

It is clear that it is no longer holomorphic in the c<j-plane. Indeed, on the lines (p = (pr^ia 
one finds a;(0) = LOy^ijl due to the quasiperiodicity property 

ci;(0 + 2m) = ci;(0) + i, uicp + 27r) = ci;(0). 

One has also 

a;r(±ia) = a;r(27r ± ia) = 0. 

The derivative ^ is double periodic with periods (27r, 2ia) and has double pole in the 
fundamental domain 

< 3?e0 < 2tt, -a < S>m0 < a. 

Therefore it must have just two zeros giving two extremal points of a;(0): one minimum 
of cUr located at (pmin + and one maximum located at 27r — (pmin + ^o. Both these 
extrema are considered with respect to line QuKf) — a, and (pmin £ (0,27r). Let us 
denote Qq = cj,.(27r — (f)min + ^a). Then it is evident that the function (j){uj) should have 
two cuts in the tu-plane represented by the segments —Qq < uJj. < f^O) = 1/2 and 
—flo < cUr < ^0, — ~l/2- Following Haldane [5], let us express as a Cauchy integral 
over the contour along the image of the boundary of the fundamental domain and use the 
symmetry properties of a;(0). We skip these rather long but in fact simple considerations. 
Only integral over the finite interval remains and after integrating by parts we obtain 

0(u;) = r dxpo{x)-ln (10) 

The still unknown function Pq{x) is normalized due to the properties of the function 0: 
a;(0 + 27r) = a;(0) and the integral representation (10), 

r° p^{x)dx = 1. (11) 
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On the other hand, we know that 0i(a; ± i/2) = ±ia for all real cu in the interval — < 
CO < D,o- This gives an integral equation for the function po{x) entering (10) of the form 
quite similar to (3), 



Note also that the same equation can be obtained with the use of quasiperiodicity property 
u!{(f) + 2ia) = u){4>) + i and the representation (10). 

The equations (3) and (12) become completely identical if one puts A = Qq, i.e. the 
meaning of the parameter A is that it defines the maximal value of u!r{4>) on the segment 
< < 27r, 0j = ia due to the uniqueness of the solution of (3) mentioned above. 

It is straightforward now to verify by differentiating (10) with respect to cu that the 
derivative % has the integral representation 



Comparing both sets of formulas (3-5) and (11-13), one can easily see that the expressions 
for the spectral density of the Lax matrix in the classical limit coincide after identification 
p{x) — po{x). This completes the proof. 
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